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, Strangeness flavor has turned out to be a very effective diagnostic tool of relativis- 

tic heavy ion physics. The absolute yield provides information about conditions 
arising in first instants of the interaction. Strange hadrons are abundant allowing 
a precise study of the chemical freeze-out conditions in the dense matter fireball. 
The thermal and chemical freeze-out universality seen in strange hadrons confirms 
C^") ' sudden hadronization as the breakup mechanism. A plausible cause of sudden fire- 

ball breakup is the mechanical instability arising when a quark-gluon plasma phase 
supercools deeply in its adiabatic expansion. Applying these ideas, we interpret 
CERN-SPS and the first results on hadronic particle ratios obtained at RHIC. 

> 

in 

£NJ , 1 Introduction 

o 

It is believed today that a new state of matter has been formed in relativistic 
nuclear collisions at CERN. The question is if this is the hot quark matter 
\ (quark-gluon plasma) state of matter, which has been for 30 years a predicted 

— «« ■ new form of elementary matter: as early as 1970, the development of the 

'"^ \ quark model lead to a first consideration of quark matter stellar structure, tl 

These ideas were deepened by the development of the quantum many body 
theory of quark matter, au which lead on to the formal recognition within the 
framework of asymptotically free quantum-chromodynamics (QCDi-that a 
very high temperature perturbative quark matter state must exist. □ Other 
reasoning based on a study of the 'boiling state' of dense hadron gas within 
the scheme of Hagedorn's statistical bootstrap of hadronic matter has inde- 
■ pendcntly lead from a different direction to the consideration of the transition 

to a hadron substructure phase. □ Considering the present day lattice-QCD 
numerical simulations, a we understand today quite well the properties of the 
baryon number free quark-gluon plasma (QGP), as we today call hot quark 
matter 

However, it is still common to study ad-hoc ideas about new phases of 
fundamental matter. This is a natural approach considering the history of 
nuclear physics which developed without fundamental knowledge about the 
quark/QCD nature of nuclear interactions. Nuclear physics advanced by the 
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development of models, and this tradition remains strong. In fact there was 
little initial interest in quarks, gluons, and QCD among nuclear physisists. 
The first of the series of formative workshops in the field of relativistic heavy 
ion collisions, the "Bear Mountain" meeting in November 1974 had not a single 
mention of quarks, let alone of quark matter. At the time ultra dense nuclei, 
multi-hyperon nuclear states, and "pion condensates", were all considered 
as the discovery potential of these new and coming tools of nuclear physics 
research. 

As the ideas about QGP formation in high energy nuclear collision ma- 
tured, an unexpected challenge emerged 20 years ago: how can the locally 
deconfined state which exists a mere 10 _22 s be distinguished from the gas of 
confined hadrons? This is also a question of principle, considering that both 
quark and hadron pictures of the reaction could be equivalent, for it has been 
argued that a quark-gluon based description is merely a change of Hilbert 
space expansion basis, if the rules of quantum mechanics with pure states are 
considered. This view contradicts the intuition, and about everybody in our 
field has come to believe that QGP phase is observable, since the quantum 
solution decoheres, even if the argument is being made with a lot of hand 
waving. 

In the Galilean tradition, such a difficult question about existence and 
observability of a new phase of elementary matter, the QGP, must be settled 
by an experiment. This requires that a probe of QGP operational on the 
collision time scale of 5T0 -23 s is developed, which is sensitive to the local color 
charge deconfinement, and that it depends on the gluon degree of freedom, 
which is the characteristic new ingredient of the quark matter phase. 

One of us (JR) proposed strangeness as the signature of QGP noticing 
more than 20 years ago that when color bonds are broken, the chemically 
equilibrated deconfined state has an unusually high number of strange quark 
pairs. A study of the dynamics of strangeness (chemical) equilibration pro- 
cess has shown that only the gluon component in the QGP is capable to 
produce strangeness rapidly Q, allowing the formation of (nearly) chemically 
equilibrated dense phase of deconfined, hot, strangeness-rich quark matter in 
relativistic nuclear collisions. Therefore, strangeness enhancement is related 
directly to the presence of gluons in QGP. 

The high density of strangeness available in the-_fireball favors during 
hadronization the formation of multi-strange hadraos, uii3 which are produced 
quite rarely if only individual hadrons collide. I 11 " 12 " A large enhancement of 
multi strange (anti)baryons has been continuously advanced and defended by 
one of us ( JR) as a characteristic signature of the QGP over the past 20 years. 
A systematic enhancement has been now reported, rising with strangeness 
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content. t2l Here the enhancement of strange antibaryons is reported with 
a base being the expectation derived from the scaled nucleon-nucleon and 
nucleon- nucleus reactions. 

These experimental results are consistent with particle production occur- 
ring from a strangeness dense source in which the strange quarks are already 
made, are freely moving around and are readily available. This is by defi- 
nition, the deconfined state we refer to as the quark-gluon plasma. We are 
not aware of a consistent interpretation of the rich and diverse experimen- 
tal results not using QGP formation. It is important here to remember that 
since there are many results on strange hadron production, any alternative 
description has to be tested on all data available. 

We think that strangeness enhancement has turned out to be a very prac- 
tical signature of quark-gluon plasma physics. This report summarizes our 
work of past 18 months in -this, area - for further details the reader may con- 
sult our recent articles. EJll3'E^lllZl In the following section || we will describe 
our most recent advances in the study of properties of QGP at time of its for- 
mation and its sudden breakup. lJ The initial conditions are determining the 
overall strangeness yield. We use the properties of QGP at sudden breakup 
in the following section ||, testing the consistency of thiixhemical freeze-out 
condition with the properties of the QGP. EJ Recently, lij we have been able 
to understand the mechanism of sudden QGP hadronization in terms of a 
mechanical instability arising within a supercooled and rapidly expanding 
fireball, as we will discuss in section ||. 

An important experimental observation pointing to sudden QGP bteakup 
is that strange baryons and antibaryons have-very similar mj_-spectra Ii3 and 
are thus produced by the same mechanism. 113 Moreover, as we will show in 
section ||, these spectra imply that the thermal freeze-out of these rarely pixu 
duced particles occurs together with the chemical (abundance) freeze-out. liH 
This occurs for all collision centralities measured by experiment WA97. 

In recent months first strangeness results have been obtained at RHIC. 
Our analysia_effprt regarding these results has just begun. Some of our RHIC 
predictions E2IO can already be compared with experiment, and we will show 
in section |^ to be in agreement with our picture of how strangeness and QGP 
evolves at RHIC. 



2 QGP fireball 
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2.1 Formulation of the QGP model 

Using the latest lattice-llCD results, we have been able to better_model the 
properties of the QGP, Ej compared to our earlier considerations. E3 The key 
ingredients of our current approach are: 



1. We relate the QCD scale to t 
scale the Matsubara frequency 
potential, it is the QCD scale) 



i temperature T = 1/(3, we use for the 
(fj, without a subscript is not a chemical 



This extension to finite chemical potential /i q , or equivalently quark fu- 
gacity A q = exp/i q /T, is motivated by the form of plasma frequency 
entering the computation of the vacuum polarization function. E-3 

We obtain the interacting strength a s (fi) integrating numerically the 
renormalization group equations 

^ = - b o* 2 s-b 1 al + ...^f3r t - (2) 

/32° rt is the beta-function of the renormalization group in two loop ap- 
proximation, and 

ll-2n f /3 51-19nf/3 
bo ~ 2^ ' h = 4^ ■ 

/?2 Crt does not depend on the renormalization scheme, and solutions of 
Eq. (|^) differ from higher order renormalization scheme dependent results 
by less than the error introduced by the experimental uncertainty in the 
measured value of a s ((i = M z ) = 0.118 + 0.001 - 0.0016. 

We introduce, in the domain of freely mobile quarks and gluons, a finite 
vacuum energy density: 

* = 0.19^X 
fm 3 

This also implies, by virtue of relativistic invariance, that there must 
be a (negative) associated pressure acting on the surface of this volume, 
aiming to reduce the size of the deconfined region. These two properties 
of the vacuum follow consistently from the vacuum partition function: 

\nZ v&c = -BV(3. (3) 
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The partition function of the quaxk-gluon liquid comprises interacting 
gluons, n q flavors of light quarks, and the vacuum i3-term. We incor- 
porate further the strange quarks by assuming that their mass in effect 
reduces their effective number n s < 1: 
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Properties of QGP-liquid 

We show properties of the quark-gluon liquid defined by Eq.|3|in a wider range 
of parameters at fixed entropy per baryon S/B, in the range S/B = 10-60 in 
step of 5 units. In the top panel in figure ^, we show baryo-chemical potential 
/it,, in middle panel baryon density u/uq, here Uq = 0.16/fm 3 , and bottom 
left the energy per baryon E/B. In top and middle panel the low entropy 
results are top-left in figure, in bottom panel bottom left. The highlighted 
curve, in figure [l], is for the value S/B = 42.5, applicable to the 158 ^4GeV 
Pb-Pb collisions. The dotted line, at the minimum of E/B\g/ B , is where the 
vacuum and quark-gluon gas pressure balance. 

Little entpapy is produced during the nearly adiabatic evolution of the 
QGP fireball. 123 Thus the lines in the lower panel of figure [l] characterize the 
approximate trajectory in time of the fireball. After initial drop in energy 
per baryon due to transfer of energy to accelerating expansion of the fireball, 
during the deep supercooling process, the motion is slowed and thus energy 
per baryon increases. The thick line is our expectation for the fireball made 
in Pb-Pb interactions at the projectile energy 158^4 GeV. The cross_shows 
the result of chemical freeze-out analysis presented in next section ^. EJ 

We also have compared the chemical freeze-out conditions with the phase 
transition properties. The hadron gas behavior is obtained evaluating and 
summing the contributions of all known hadronic resonances considered to 
be point particles. The thin solid line in the T, /ib plane in figure |^ shows 
where the pressure of the quark-gluon liquid equals the equilibrated hadron 
gas pressure. When we allow for finite volume of hadrons, E3 we find that 



tab20002: submitted to World Scientific on February 8, 2008 5 



Q Lj I I I I I I I I I I I I I I I I I L 

0.1 0.2 0.3 

T [GeV] 

Figure 1. From top to bottom: fjij,, n/no and E/B; lines shown correspond to fixed 
entropy per baryon S/B = 10 to 60 by step of 5 (left to right). Thick solid lines: result 
for S/B = 42.5. Limits: energy density e qi9 = 0.5GeV/fm 3 and baryo-chemicalrpptential 
ptf, = 1 GeV. The experimental points denote chemical freeze-out analysis result. 



the hadron pressure is slightly reduced, leading to some (5MeV) reduction 
in the equilibrium transition temperature, as is shown by the dashed line in 
figure ^. For vanishing baryo-chemical potential, we note in figure ^ that 
the equilibrium phase transition temperature is T pt ~ 172 MeV, and when 
finite hadron size is allowed, Tf p ~ 166 MeV, The scale in temperature we 
discuss is result of comparison with lattice gauge results. Within the lattice 
calculations it arises from the comparison with the string tension. 
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Figure 2. Thin solid and dashed lines: equilibrium phase transition from hadron gas to 
QGP liquid without and with excluded volume correction, respectively. Dotted: breakup 
condition at shape parameter k = 0.6, for expansion velocity = 0, 1/10, 1/6, 1/5, 1/4 and 
1/3, and thick line for v c =p-0.54. The experimental point denotes chemical non-equilibrium 
freeze-out analysis result. Ej 



The dotted lines, in figure |[ show where the flowing QGP liquid has 
a pressure which balances the vacuum pressure, and thus these lines cor- 
respond to sudden break up at a velocity for (from right to left) v% = 
0, 1/10, 1/6, 1/5, 1/4 and 1/3, see section |], specifically the dotted lines in fig- 
ure p| co rrespond to the condition Eq. (|^) using the shape parameter k = 0.6, 
Eq. (|12|). The last dotted line corresponds thus to an expansion flow with 
the velocity of sound of relativistic noninteracting massless gas. The thick 
solid line corresponds to an expansion with v c — 0Ji4. The hadron anal- 
ysis result is also shown as presented in section || EJ Comparing in figure 
U thin solid/dashed with the thick line, we recognize the deep supercooling 
as required for the explosive fireball disintegration. The super-cooled zero 
pressure P = baryonfree QGP temperature is at T sc = 157 MeV, (see the 
intercept of the first dashed line to the right in figure ^) and an expanding 
fireball can deeply super-cool to Td sc — 147 MeV (see the intercept of thick 
solid line) before the mechanical instability occurs. 
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2.3 Initial conditions in QGP formation 

The formation of the QGP occurs well before chemical equilibration of quarks. 
We model this situation varying rif. For nf = 1 we show, in figure |[ lines 
of fixed energy per baryon E/B =3, 4, 5, 6, 8, 10, 20, 50 and 100 GeV, The hor- 
izontal solid line is where the equilibrated hadronic gas phase has the same 
pressure as QGP-liquid with semi-equilibrated quark abundance. The free en- 
ergy of the QGP liquid must be lower (pressure higher) in order for hadrons 
to dissolve into the plasma phase. The dotted lines in figure ||, from bottom 
to top, show where the pressure of the semi-equilibrated QGP phase is equal 
to n = 20%, 40%, 60%, 80% and-,100%, i] being the 'stopping' fraction of the 
dynamical collisional pressure c3: 

P col = VPo^- (6) 

The rationale to study, in figure | lines at fixed E/B is that, during 
the nuclear collision which lasts about 2Rn /jl"2c ~ 13/18fm/c, where jl 
is the Lorentz factor between the lab and CM frame and Rn is the nuclear 
radius, parton collisions lead to a partial (assumed here to be 1/2) chemi- 
cal equilibration of the hadron matter. At that time, the pressure exercised 
corresponds to collisional pressure P co i . This stopping fraction, seen in the 
transverse energy produced, is about 40% for S-S collisions at 200A GeV and 
60% for Pb-Pb collisions at 158 A GeV. If the momentum-energy and baryon 
number stopping are similar, as we see in the experimental data, then the SPS 
collisions at 160-200AGeV are found in the highlighted area left of center of 
the figure. In the middle of upper boundary of this area, we would expect 
the beginning evolution of the thermal but not yet chemically equilibrated 
Pb-Pb fireball, and in the lower left corner of the S-S fireball. We note that 
the temperature reached in S-S case is seen to be about 25 MeV lower than in 
the Pb-Pb case. The lowest dotted line (20% stopping) nearly coincides with 
the non-equilibrium phase boundary (solid horizontal line, in figure [|) and 
thus we conclude that this is, for the condition rif = 1, the lowest stopping 
that can lead to formation of a deconfined QGP phase. Such a low stopping 
would be encountered possibly in lighter than S-S collision systems or/and at 
large impact parameter interactions of larger nuclei. 

The highlighted area, right of center of the figure ||, corresponds to the 
expected conditions in Pb-Pb collisions at AQA GeV. If we assume that the 
stopping here is near 80%, then the initial conditions for fireball evolution 
would be found towards the upper right corner of this highlighted area. We 
recognize that the higher stopping nearly completely compensates the effect 
of reduced available energy in the collision and indeed, we expect that we form 
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Figure 3. Contours of energy per baryon in QGP in the T—\ q plane for nf = 1: From right 
to left E/B=3, 4, 5, 6, 8, 10, 20, 50 and 100 GeV . Thin, nearly horizontal line: hadronic 
gas phase has the same pressure as the QGP-liquid with semi-equilibrated quark flavor. 
Dotted lines from bottom to top: pressure in QGP liquid equals 20%, 40%, 60%, 80%, and 
100%, of the dynamical collisional pressure. 



QGP also at these collision energies. It is important to realize that we are 
entering a domain of parameters, in particular X q , for which the extrapolation 
of the lattice results is not necessarily reliable, and thus our equations of state 
have increased systematic uncertainty. 

3 Fermi-model data analysis 

3. 1 Properties of hadronic matter fireball 

A full account of. our prior analysis of the 158^4 GeV Pb-Pb collision system 
has appeared, cj We briefly summarize the results that we require for the 
study of QGP properties at fireball breakup. In table [j], in upper section, 
we present the parameters Tf (the chemical freeze-out temperature), v c (the 
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collective flow velocity at sudden breakup), X q (the light quark fugacity), A s 
(the strange quark fugacity), "f q (the light quark phase space occupancy), 7 S 
(the strange quark phase space occupancy). These chemical properties are 
derived from analysis of all hadrons excluding fl and fi, which data points 
are not following the same systematic production pattern. These parameters 
characterize completely the physical properties of the produced hadrons, and 
these properties are shown in the bottom section of table [j]. 

In the heading of the table, the total error, x 2 is shown, along with 
the number of data points N, parameters p and data point constraints r. 
The confidence level that is reached in our description is near or above 90%, 
depending on scenario considered. The scenarios we consider are seen in the 
columns of table |l|: an unconstrained description of all data in the first column, 
constraint to exact strangeness conservation in the observed hadrons, second 
column. Since in both cases the parameter "f q assumes value that maximizes 
the entropy and energy content in the pion gas, we assume this value in the 
so constrained third column. 

We can now check the consistency between the statistical parameters (top 
panel of table [l]) and the physical properties of the fireball (bottom panel of 
table |l|) which are maintained in the process of hadronization. We note that 
the energy shown in this table, is the intrinsic energy in the flowing frame. 
The CM-laboratory energy includes the kinetic energy of the flow and thus is 
greater, to be obtained multiplying the result shown in table [| by the Lorenz 
factor 7=1/ \/i — u 2 = 1.19. Thus the initial value of the energy per baryon 
that the system had before expansion started has been E°/B ~ 9.3 GeV. 

3.2 Interpretation in terms of QGP 

In the bottom panel in figure [j], we saw that the Temperature Tf = 143 ± 
3MeV and intrinsic (frame at rest) energy per baryon E/B = 7.8 GeV where 
just at S/B — 42.5 seen tabic |l} Similarly, in the top panel, the baryo-chemical 
potential /if, = 3Tf\n\ q — 204 ± lOMeV is as required for the consistency of 
QGP properties. The fireball hadron freeze-out and/or breakup condition is 
thus found well below the QGP to HG phase transition temperature. Both 
the specific energy and entropy content of the fireball are consistent with the 
statistical parameters Tf and /if, according to our equations of state of the 
quark-gluon liquid. We can also evaluate the hadronic phase space energy 
density. We allow the excluded volume correction. E3 Considering that the 
point particle phase space energy density e p t = 1.1 GcV/fm 3 , we obtain shg — 
0.4GeV/fm 3 at chemical freeze-out, using the value of B = 0.19GeV/fm 3 . 
This energy density is neraly the same as one finds in QGP. Thus we see that 
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Table 1. Results of study of Pb— Pb hadron production: in the heading: the total 
quadratic relative error x^i number of data points N, parameters p and redundancies r; in 
the upper section: statistical model parameters which best describe the experimental results 
for Pb-Pb data. Bottom section: specific energy, entropy, anti-strangeness, net strangeness 
of the full hadron phase space characterized by these statistical parameters. In column one, 
all statistical parameters and the flow vary. In column two, we fix \ s by requirement of 
strangeness conservation, and in column three, we fix -f q at the pion condensation point 

7q = lq- 





Pb|„ 


Pb|f 


Pb|- 


Xt5 N;p;r 


2.5; 12;6;2 


3.2; 12; 5; 2 


2.6; 12; 5; 2 


T f [MeV] 


142 ± 3 


144 ± 2 


142 ± 2 


v c 


0.54 ± 0.04 


0.54 ± 0.025 


0.54 ± 0.025 


\ 


1.61 ± 0.02 


1.605 ± 0.025 


1.615 ± 0.025 


A s 


1.09 ± 0.02 


1.10* 


1.09 ± 0.02 


lq 


1.7 ± 0.5 


1.8± 0.2 


7 c* _ e m„/2T f 


Is/lq 


0.79 ± 0.05 


0.80 ± 0.05 


0.79 ± 0.05 


Ef/B 


7.8 ± 0.5 


7.7 ± 0.5 


7.8 ± 0.5 


Sf/B 


42 ± 3 


41 ± 3 


43 ± 3 


Sf/B 


0.69 ± 0.04 


0.67 ± 0.05 


0.70 ± 0.05 


0/ - s f )/B 


0.03 ± 0.04 


0* 


0.04 ± 0.05 



the chemical freeze-out is due to the direct fireball breakup. 

More generally, the chemical analysis evidence for the QGP nature of the 
fireball is as follows: li-3 

a) The value of strange quark fugacity A s we find in our analysis can also 
be obtained from the requirement that strangeness balances, (N s — N s ) = , 
in QGP, which for a source in which all s an s quarks are unbound and 
thus have symmetric phase space, implies X s — 1 . However, the Coulomb 
distortion of the strange quark phase space plays an important role in the un- 
derstanding of this constraint for Pb-Pb collisions, citeRaf91 leading to the 
Coulomb-deformed value A s ~ 1.1, which is identical to the value obtained 
from experimental data analysis A" = 1.09 ± 0.02 . 

b) The phase space occupancy of light quarks 7 g even if it were before gluon 
fragmentation near or at the equilibrium value j q = 1, must increase above 
this equilibrium value. There is an upper limit: 7 9 < 7^ = e m "^ 2Tf ~ 1.67, 
which arises to maximize the entropy density in the confined hadron phase. 
In our analysis we find values in this narrow interval 1 < 7 g < 7^ 

c) The strange quark phase space occupancy 7 S we compute within the frame- 
work of kinetic theory where strangeness pair production arises in gluon fu- 



tab20002: submitted to World Scientific on February 8, 2008 11 



sion, B. The result depends on temperature reached in early stages of the col- 
lision, and following dilution effect in which the already produced strangeness 
can even over saturate the 'thiner' low temperature phase space. Moreover, 
some gluon fragmentation also enriches 7 S as measured by hadron abundance, 
and thus we find gmma s > 1. We note that the ratio j s /jq has been earlier 
confounded with 7 S , which therefore is stated in table |l|. 
d) The strangeness yield, N s /B = N s /B ~ 0.7, predicted early on as the re- 
sult of QGP formation □ is also one of the results of data analysis seen in table 
[I]. This result we_obtain within the modern kinetic model of strangeness pro- 
duction in QGP. lia The results are shown in figure |] as function of the intrinsic 
specific energy available in the fireball E/B, for the three collision systems S- 
Au/W/Pb (short-dashed line), Ag-Ag (long dashed) and Pb-Pb (solid line). 
The solid square is the 'experimental' result obtained in our analysis of the 
S-Au/W/Pb system, and the open square is for Pb-Pb. The experimental 
value of stopping we employ (see legend to figure ^) were obtained by NA49 
collaboration, tB. 

Stopping enters the computation of the initial temperature from the colli- 
sion energy, see Eq. (@). We extrapolate as function of E/B, assuming that the 
stopping fraction for the collisional pressure is known, and as here assumed, 
not changing with collision energy in the narrow interval shown in figure ^, 
leading to a steep rise of strangeness production with collision energy. 

4 Supercooling and sudden breakup 

4-1 Instability condition 

Sudden hadronization of dense matter fireballs formed in 158^4 GeV Pb-Pb 
collisions requires that a global instability of the dense hadronic matter fire- 
ball is encountered in its evolution. Given the hadronization temperature 
we found, T ~ 143 MeV, which is about 15MeV below the equilibrium phase 
transition of QGP phase to hadron gas it is a natural step to consider the pos- 
sibility of supercooling of the QGP phase. As many of us have experienced in 
diverse daily life experiences, once the supercooled phase is disturbed, a truly 
sudden break-up occurs. QGP instabiliti^was suggested by Csorgo and Cser- 
nai Ea to explain sudden hadronizatiori £3 Despite extensive ensuing study of 
potential microscopic mechanisms, l 31 H 3 4P 3 l a convincing elementary instability 
mechanism has not been found. 

As the fireball expands, energy is transfered from thermal motion to col- 
lective flow in adiabatic fashion. At a certain point in deep supercooling, 
energy flow is reversed, and the matter flow begins to slow down. However, if 
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Figure 4. QGP fireball specific per baryon strangeness abundance as function of E/B 
energy per baryon in the fireball, for running a s (Mz) = 0.118 and m B (lGeV) = 200 MeV, 
corresponding to m B (Mz) = 90 MeV. Solid lines: for Pb— Pb, stopping 57% and 360 par- 
ticipants, long-dashed lines: for Ag-Ag, stopping 54% and 180 participants, short-dashed 
lines: for S-Au/W/Pb stopping 52% and 90 participants in QGP fireball. The solid square 
is the result of an analysis for S-Au/W/Pb 200A GeV reaction system, and open square for 
Pb-Pb at 158A GeV, as shown in the lower section of table ttL 



the flow velocity remains large but the surface pressure V — > 0, continued ex- 
pansion must rip the fireball apart. This is the natural mechanical instability 
we proposed as the sudden hadronization condition. Ii3 This simple mechanism 
does not require a microscopic process for the onset of instability. 

We consider the exploding fireball dynamics in its center of momentum 
frame of reference. The surface normal vector of exploding fireball is ft, and 
the local velocity of matter flow v c . The rate of momentum flow vector V at 
the surface is obtained from the energy-stress tensor T^V- Ej 

p = pWf?+(pW+ e W)^M. (7) 

1 — 1 ) * 
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The upper index (i) refers for the intrinsic energy density e and pressure 
P of matter in the frame of reference, locally at rest, i.e., observed by a co- 
moving observer. We omit the superscript (i) in the following. For the fireball 
expansion to continue, V = \V\ > is required. For V — > at v c =/= 0, we 
have a conflict between the desire of the motion to stop or even reverse, and 
the continued inertial expansion. 

When the flow velocity remains large but V — > 0, the intrinsic pressure P 
must be negative. As illustration consider the fireball to be made of a quark- 
gluon liquid confined by an external vacuum pressure B. The total pressure 
and energy comprise particle (subscript p) and the vacuum properties: 

P = P P -B, e = e p + B. (8) 

Eq. (0) with V = thus reads: 

V V 'ft 

Bn = P p n+(P p + e p )-^ Y , (9) 

and it describes the (equilibrium) condition where the pressure of the expand- 
ing quark-gluon fluid is just balanced by the external vacuum pressure. 

Expansion beyond V — > is in general not possible. A surface region 
of the fireball that reached it but continues to flow outwards must be torn 
apart. This is a collective instability and thus the ensuing disintegration 
of the fireball matter will be very rapid, provided that much of the surface 
reaches this condition. We adopt the condition V = at any surface region 
to be the instability condition of an expanding hadron matter fireball. 

Negative internal pressure P < is a requirement. At this stage the fire- 
ball must thus be significantly supercooled. The adiabatic transfer of internal 
heat into accelerating flow of matter provides the mechanism which leads on 
the scale of t = 2 10~ 23 s to the development of this 'deep' supercooling. 

4-2 Consistency with QGP properties 

The sudden break up due to mechanical instability of course will occur for 
all phases of elementary matter and thus we consider next in how far the 
observed freeze-out conditions and other features of hadronization process 
are consistent with the QGP mechanical breakup, and explore if experiments 
show evidence for supercooling. 

It is possible to determine experimentally if the condition P < has been 
reached. Namely, the Gibbs-Duham relation for a unit volume: 

P = Ta + n b v b -e, (10) 



tab20002: submitted to World Scientific on February 8, 2008 14 



relates the pressure, to entropy density a = S/V , energy density e = E/V , 
and baryon density v\> — b/V , V is the volume, T is the temperature, and /ib 
the baryochemical potential. Dividing by e we obtain: 

£ E/S E/B 1 J 

The microscopic processes governing the fireball breakup determine how the 
quantities entering the right hand side of Eq. ( pi] ) are changed as hadrons 
emerge. Understanding this we can determine, if the intrinsic fireball pressure 
prior to breakup, has been negative. 

The energy E and baryon content b of the fireball are conserved. Entropy 
S is conserved when the gluon content of a QGP fireball is transformed into 
quark pairs in the entropy conserving process G + G —> q + q. Similarly, when 
quarks and antiquarks recombine into hadrons, entropy is conserved in the 
range of parameters of interest here. Thus also E/B and S/B is conserved 
across hadronization condition. The sudden hadronization process also main- 
tains the temperature T and baryo-chemical potential /Zh across the phase 
boundary. What changes in QGP breakup are the chemical occupancy pa- 
rameters. As gluons convert into quark pairs and hadrons j g — > but the 
number occupancy of light valance quark pairs increases 7 g > 7 go ~ 1 in- 
creases significantly, along with the number occupancy of strange quark pairs 

Is > 7s ^ 1. 

Evaluating Eq. ( pi] ) using the results of the data analysis presented above, 
we indeed obtain Pf < 0. The magnitude of |P/| can vary between a few per- 
cent (in_lerms of energy density Ej/V), up to 20% for the latest published 
result. M The precise value, which arises from several cancellations of larger 
numbers is sensitive to the strategy of how the currently available experimen- 
tal data is described, e.g., if strangeness conservation is implemented, and if 
so, if differentially at each rapidity, or as an overall conservation law; how 
many high mass resonances can be excited in hadronization process, etc ... 

Understanding in detail the breakup condition V — > requires that we 
model the shape and direction of flow in the late stage of fireball evolution, 
obviously not an easy task. However, considering n ■ V —>■ 0, we find the 
constraint: 

-PV/E =K ^L^ K =(v c -n) 2 /v 2 c . (12) 



1 + PV/E 1 - v, 

For an exactly symmetrical, spherical expansion the two vectors v c and n are 
everywhere parallel, thus n — > 1. However, in 158^4 GeV Pb-Pb reactions the 
longitudinal flow is considerably greater than the transverse flow, c£l and we 
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note k — > for a longitudinally evolving cylindrical fireball. For the Pb-Pb 
collisions considered here, our analysis suggest 0.1 < k < 0.6. 

We now substitute, in Eq. ( |l2|) , the fireball matter properties employing 
the Gibbs-Duham relation, Eq. (|10|), and arrive at: 

§-(^){ 1 + T^}- M 

Eq. (|l^) establishes a general constraint characterizing the fireball breakup 
condition. 

Deep supercooling requires a first order phase transition, and this in turn 
implies presence of latent heat B. Physical consistency then requires presence 
of external (negative) vacuum pressure —B. We now combine the theoretical 
properties of the QGP equations of state with the dynamical fireball properties 
in order to constrain B. Reviewing Eq. (11), we obtain: 

PV 

-—e QGP + P p =B, (14) 

To evaluate B, we note that lattice results for £qgp are well represented by 
£qgp = aT 4 , with a ~ 11, value extrapolated for_the number of light quark 
flavors being rif = 2.5 at the hadronization point. We obtain, for the fireball 
formed in Pb-Pb reactions, 

^4 



0.2 • 117£ ~ 0.17 GeV/fm < B. 



5 Thermal freeze-out 

5.1 Direct and decay m_i_- spectra 

We next study the shape of hadron m^-spectra. These are strongly influenced 
by resonance decays, which have also contributed in an important fashion to 
the yield of all stable hadrons observed. 

The final particle distribution is composed of directly produced particles 
and decay products of heavier hadronic resonances: 

dN x _ dN x 
dm± dm± 

Here, R(M, Mr, Y) — ► X(m, mr, y) + 2(m2) + • • •, where we indicate by the 
arguments that only for the decay particle X we keep the information about 
the shape of the momentum spectrum. 



+ dN 

direct w n .vim ^TYl 



X 



(15) 
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In detail, the decay contribution to yield of X is: 

dN x g r b [ Y + Jv f M - + 2 d 2 N R 

dY / dM T J — (16) 



dm\dy Airp* J Y J Mt _ 1 dM 2 dY 

^PJpJ - {ME* - M T m T cosh AF} 2 

We have used AF = Y — y, and ^fs is the combined invariant mass of the 
decay products other than particle X and E* = [M 2 — m 2 - m 2 .)/2M, p* = 
\JE* 2 — m 2 are the energy and momentum of the decay particle X in the rest 
frame of its parent. The limits on the integration are the maximum values 
accessible to the decay product X: 

Y± -yisinh" 1 ( -H- 
V m T 



E*m T cosh AY ±pt\/p* 2 - m% sinh 2 AY 

M T± = M ^ 

mip sinh Ay + m 2 

The theoretical primary particle spectra (both those directly produced 
and parents of decay products) are derived from the Boltzmann distribution by 
Lorenz-transforming from a flowing intrinsic fluid element to the CM-framc, 
and integrating over allowed angles between particle direction and local flow. 

We introduce in the current analysis two velocities: a local flow velocity 
v of fireball matter where from particles emerge, and hadronization surface 
(breakup) velocity which we refer to as vj 1 = dtf/dxf. Particle production 
is controlled by the effective volume element, which comprises this quantity. 
In detail: 

,n u , (, v^-A d 3 xd 3 p 

dS»p» =du,ll- 1 , du EE . (17) 

The Boltzmann distribution we adapt has thus the form 
d 2 N f_ v^-p 



dmrdy \ E 



Vf 1 • P \ _ Ett v-p\ 

1 J — — 7 m T cosh ye 7t( - e ' , (18) 



where 7 = 1/Vl — v 2 . 

The normalization for each hadron type h = X, R is 

n 

N h = V QG pI[\ i1i . 

ieh 
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We use the chemical parameters Xi,jii = q,s which are as defined in EJ and 
commonly used to characterize relative and absolute abundances of light and 
strange quarks. 

The best thermal and chemical parameters which minimize the total rel- 
ative error \t, 



for all experimental measurement points Fi having measurement error AFi are 
determined by considering simultaneously all strange baryon and antibaryon 
results of experiment WA97, lla i.e., A, A, 5, S, £1 + £1, as well as Kg, except 
that we assign additional systematic error to these data (see below). 

5.2 A-, A-, H-, S-mj_- spectra 

In recent months experiment WA97 determined the absolute normalization of 
the published m± distribution, 113 and we took the opportunity to perform the 
shape analysis together with yield analysis in order to check if the thermal 
and chemical freeze-out conditions are the same. E3 Our analysis addresses 
the data differentiated in four different collision centralities. This allows to 
study the effect of the changing volume within the limits of precision set by 
the current data. 

Since particle spectra we consider have a good relative normalization, only 
one parameter is required for each centrality in order to describe the absolute 
normalization of all six hadron spectra. This is for two reasons important: 

a) we can check if the volume from which strange hadrons are emitted grows 
with centrality of the collision as we expect; 

b) we can determine which region in mj_ produces the excess of f2 noted in 
the chemical fit. liJ 

However, since the normalization Vqgp common for all particles at given 
centrality comprises additional experimental acceptance normalization which 
we have not yet studied, we have not normalized the value of the fireball 
emission volume at each centrality. Hence we will be presenting the volume 
parameter as function of centrality in arbitrary units. 

If thermal and chemical freeze-outs are identical, our present results must 
be consistent with earlier chemical analysis of hadron yields. Since the ex- 
perimental data we here study is dominated by the shape of mj_-spectra and 
not by relative particle yields, our analysis is de facto comparing thermal and 
chemical freeze-outs. 
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Figure 5. Thermal freeze-out temperature T for different centrality bins compared to the 
chemical freeze-out analysis. 



We next show the parameters determining the shape of the m± distribu- 
tions, that is T, v, Vf. As function of the centrality bin 1,2,3,4 with the most 
central bin being 4 we show in figure [5] the freeze-out temperature T of the m± 
spectra. The horizontal lines delineate range of result of the chemical analy- 
sis. It is reassuring that we find a result consistent path the purely chemical 
analysis of data that included non-strange hadrons. LL-3 There is no indication 
of a significant or systematic change of T with centrality though there is a 
somewhat higher value in the most central 4th bin. This is consistent with the 
believe that the formation of the new state of matter at CERN is occurring in 
all centrality bins explored by the experiment WA97. Only most peripheral 
interactions produce a change in the pattern of strange hadron production. 

The (unweighted) average of all results shown in figure produces a freeze- 
out temperature at the upper boundary of the the pure chemical freeze-out 
analysis result, T ~ 145 MeV. It should be noted that in chemical analysis we 
have not distinguished between the flow v and freeze-out Vf velocity and both 
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Figure 6. Thermal freeze-out flow velocity v (top) and break up (hadronization) velocity 
Vf for different centrality bins. Upper limit Ve = 1 (dashed line) and chemical freeze-out 
analysis limits for v (solid lines) are also shown. 



are denoted as v c , which may be the cause of this slight difference between 
current analysis average and the earlier purely chemical analysis result. 

The magnitudes of the collective expansion velocity v and the break-up 
(hadronization) speed parameter Vf are presented in figure ^. For v (lower 
part of the figure) we again see consistency with earlier chemical freeze-out 
analysis results, and there is no confirmed systematic trend in the behavior 
of this parameter as function of centrality. 

Though within the experimental error, one could argue inspecting figure 
^ that there is systematic increase in transverse flow velocity v with centrality 
and thus size of the system. This is expected, since the more central events 
comprise greater volume of matter, which allows more time for development 
of the flow. Interestingly, it is in v and not T that we find the slight change 
of spectral slopes noted in the presentation of the experimental data. £3 

The value of the beak-up (hadronization) speed parameter Vf shown in 
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Figure 7. Hadronization volume (arbitrary units) for different centrality bins. 



the top portion of figure |6| is near to velocity of light which is consistent with 
the picture of a sudden breakup of the fireball. This hadronization surface 
velocity Vf was in the earlier chemical fit set to be equal to v, as there was 
not enough sensitivity in purely chemical fit to determine the value of vj. 

Unlike the temperature and two velocities, the overall normalization of 
hadron yields, V h is found to be strongly centrality dependent, as is seen in 
figure ^. This confirms in quantitative way the believe that the entire available 
fireball volume is available for hadron production. The strong increase in the 
volume by factor six is qualitatively consistent with a geometric interpretation 
of the collision centrality effect. Not shown is the error propagating from the 
experimental data which is strongly correlated to the chemical parameters 
discussed next. This systematic uncertainty is another reason we do not 
attempt an absolute unit volume normalization. 

The four chemical parameters X q , X s , 7 g , ^s/lq are shown in the following 
figures ||, ||. These parameters determine along with V h the final particle 
yield. Since we have 5 parameters determining normalization of 6 strange 
hadron spectra, and as discussed we reduce the statistical weight of Kaons, 
there is obviously a lot of correlation between these 4 quantities, and thus the 
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Figure 8. Thermal analysis chemical quark fugacity \ q (top) and strange quark fugacity X s 
(bottom) for different ccntrality bins compared to the chemical freeze-out analysis results. 



error bar which reflects this correlation, is significant. 

The chemical fugacities \ q and \ s shown in figure ||| do not exhibit a sys- 
tematic centrality dependence. This is consistent with the result we found 
for T in that the freeze-out properties of the fireball are seen to be for the 
temperature and chemical potential values independent of the size of the fire- 
ball. Comparing to the earlier chemical freeze-out result in figure || one may 
argue that there is a systematic downward deviation in X q . However, this 
could be caused by the fact that the_chemical freeze-out analysis allowed for 
isospin-asymmetric S _ (<fss) yield, Ej while out present analysis is not yet 
distinguishing light quarks. 

The ratio 7 s /7 9 shown in bottom portion of figure |^ is systematically 
smaller than unity, consistent with many years of prior analysis: when jg = 1 
is tacitly chosen, this ratio is the value of 7 S in analysis of strange baryons. 
We have not imposed a constrain on the range of "f q (top of figure and thus 
values greater than the pion condensation point 7* = e m ^/ 2T ~ \ 55 (thick 
line) can be expected, but in fact do not arise. 
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Figure 9. Thermal analysis chemical quark abundance parameter 7 9 (top) and "fs/'fq 
(bottom) for different centrality bins compared to the chemical freeze-out analysis. Thick 
line: upper limit arising with pion condensation. 



It is important to explicitly check how well the particle m^-spectra are 
reproduced. We group all bins in one figure and show in figures [h], [ll], |l2|, |l3| 
in sequence A, A, S, S. Overall, the description of the shape of the spectra is 
very satisfactory. 



5.3 Omega and Meson spectra 

In figure [l4| all four centrality bins for the sum Q + Q are shown. We see 
that we systematically under predict the two lowest m± data points. Some 
deviation at high m± may be attributable to acceptance uncertainties, also 
seen in the the S result presented earlier in figure [l^. We recall that there 
is a disagreement between theoritical expectation and the Omega yields. In 
the here presented thermal analysis, we see that this disagreement is arising 
at low momentum. We could include all experimental results in the present 
analysis since the omega data have a relative low statistical evidence and thus 
other hyperon data determine the behavior of il-spectrum. 

The 1ow-to_l anomaly also explains why the inverse m± slopes for f2, 
are smaller than the values seen in all other strange (anti)hyperons. One can 
presently only speculate about the processes which contribute to this anomaly. 
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Figure 10. Thermal analysis my spectra: A 
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Figure 11. Thermal analysis my spectra: A. 
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Figure 12. Thermal analysis my spectra: S 
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Figure 13. Thermal analysis my spectra: S 
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Figure 14. Thermal analysis mx spectra: f2 + Q. 



We note that the l-2s.d. deviations in the low mj_-bins of the + 17 spectrum 
translates into 3s. d. deviations from the prediction of the chemical analysis. 
This is mainly a consequence of the fact that after summing over centrality 
and m±, the statistical error which dominates Q+f2 spectra becomes relatively 
small, and as can be seen the low m± excess practically doubles the yield. 

For kaons K° the statistical errors are very small, and we find in a more 
in depth statistical analysis that they must be smaller than the systematic 
errors not considered. This is done checking the influence of the deletion of 
individual measurement points in the spectrum on the fit result, which influ- 
ences the outcome in case of kaons only. We assign in a somewhat arbitrary 
manner a systematic error that we arbitrarily choose to be 5 times greater 
than statistical error. 

With this modification we included the kaon data in the analysis which 
included the above presented strange baryon and antibaryon results. Because 
of the large error we assigned, there is relatively little influence the kaon 
spectra have on the baryon/ antibaryon analysis. In fact the kaon results show 
how the hyperon/antihyperon experimental results predict Kg-mj_ spectra. 
We describe the kaon data well, and especially so in the m± range which is 
the same as that for hyperons considered earlier, as is seen in figure |l5|. 

We now turn to the pion spectra measured by the WA98 collaboration. E3 
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Figure 15. Thermal analysis my spectra: K B . 



This is a very difficult matter, since at low momenta there are many mecha- 
nisms of pion production which we have not analyzed. Moreover, we have not 
fully understood how the longitudinal flow influences through the secondary 
decays the observed pion spectra. Compared to strange hadrons we study the 
7T° data in a very wide range of m± in which the yield changes by 6 orders of 
magnitude. 

We use a set of parameters obtained fro the hyperon data, i.e., Tt = 
143 MeV, setting also freeze-out surface velocity 0.99c. We include in the 
results shown in figure |l6| aside of directly produced pions the two body decay 
of the p. We allow a direct hard parton QCD component contribution, E3 of 
the form 

Ed 3 N/d 3 poc 1/pl. 

Thus, we vary four parameters: two normalizations (thermal and QCD direct 
component), and also K = 5.6 ± 1.2 , v c = 0.55 ± 0.02, with resulting best 
fits as stated. This procedure yields a good description of the data as seen 
in figure |l6|, with % 2 /dof < 1.4. Considering that only statistical error is 
being considered, and we were not evaluating contributions of 3-body decay 
resonances, we see this as a noteworthy confirmation that the hard pions are 
also created primarily in the sudden hadronization process. We note that the 
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Figure 16. Central rapidity data of experiment WA98 for n° are compared to the spectra 
expected in a sudden freeze-out reaction picture. See text for details. 



direct QCD component is at 1-20% level, and thus before proceeding further 
a more precise model of this contribution needs to be_developed. 

Our analysis is not contradicting results shown in E3, for these authors did 
not consider that the freeze-out surface velocity is different from the flow ve- 
locity, and they did not allow for direct parton-parton scattering contributions 
in their analysis of pion spectra. 

Our results confirm the finding that the apparent temperature hierarchy 
for different mass particles is due to a collective expansion of the source. EJ 

6 Strangeness at RHIC 

6. 1 Expectations 

The abundances of particles produced from QGP within the sudden freeze-out 
model are controlled by several parameters we restate here: the light quark 
fugacity 1 < X q < 1.1 , value is limited by the expected small ratio between 
baryons and mesons (baryon-poor plasma) when the energy per baryon is 
above 100 GeV, strangeness fugacity A s ~ 1 which value for locally neutral 
plasma assures that (s — s) — 0; the light quark phase space occupancy j q ~ 
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Table 2. For 7 S = 1.25, A s = 1 and y q , \ q as shown: Top portion: strangeness per baryon 
s/B, energy per baryon E / B[GeV] and entropy per baryon S/B. Bottom portion: sample 
of hadron ratios expected at RHIC. 



7<? 


1.25 


1.5 


1.5 


1.5 


1.60 


\ 


1.03 


1.025 


1.03 


1.035 


1.03 


E/B[GeV] 


117 


133 


111 


95 


110 


s/B 


18 


16 


13 


12 


12 


S/B 


630 


698 


583 


501 


571 


p/P 


1.19 


1.15 


1.19 


1.22 


1.19 


A/p 


1.74 


1.47 


1.47 


1.45 


1.35 


A/p 


1.85 


1.54 


1.55 


1.55 


1.44 


A/A 


0.89 


0.91 


0.89 


0.87 


0.89 


B-/A 


0.19 


0.16 


0.16 


0.16 


0.15 


S-/A 


0.20 


0.17 


0.17 


0.17 


0.16 




0.94 


0.95 


0.94 


0.93 


0.94 


n/s~ 


0.147 


0.123 


0.122 


0.122 


0.115 




0.156 


0.130 


0.130 


0.131 


0.122 


n/n 


1 


1. 


1. 


1. 


1. 


n+n 


0.15 


0.13 


0.13 


0.13 


0.12 


A+A 


0.19 


0.16 


0.16 


0.16 


0.15 


K+/K- 


1.05 


1.04 


1.05 


1.06 


1.05 



1.5, overabundance value due to gluon fragmentation. Given these narrow 
ranges of chemical parameters and the sudden freeze-out temperature Tf = 
150 MeV expected for nearly baryon free domain, we compute the particle 
production at break-up. 

Taking y q — I.SJ^b > we choose the value of A g , see the header of table 
^|, for which the energy per baryon (E / B) is similar to the collision condition 
(100 GeV), which leads to the range X q = 1.03 ± 0.005. We evaluate for 
these examples aside of E/B, the strangeness per baryon s/B and entropy 
per baryon S/ B as shown in the top section of the table |^. We do not enforce 
(s — s) — at central rapidity exactly, but since baryon asymmetry is small, 
strangeness is balanced to better than 2% in the parameter range considered. 
In the bottom portion of table ||, we present the compatible particle abundance 
ratios, jComputed according to the procedure we developed, and described 
above. Ilil We have given, aside of the baryon and antibaryon relative yields, 
also the relative kaon yield, which is also well determined within this approach. 

The meaning of these results can be better appreciated when we assume 
in an example the central rapidity density of direct protons is dp/dy\ CCIL t. = 25. 
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Table 3. dN/dy\ ccn t, assuming in this example dp/dy\ ccn t = 25 . 
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^0 
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n=n 


1.25 


1.03 


17 


25* 


21 


44 


39 


31 


27 


17 


16 


1.2 


1.5 
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13 


25* 


22 


36 


33 


26 


23 


13 


11 


0.7 


1.5 
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21 
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26 


23 


12 


11 


0.7 


1.5 


1.035 


18 


25* 


21 


36 


32 


26 


22 


11 


10 


0.7 


1.60 


1.03 


15 


25* 


21 


34 


30 


24 


21 


10 


9.6 


0.6 



In table |3|, we present the resulting (anti)baryon abundances. The net baryon 
density db/dy ~ 16 ± 3, there is baryon number transparency. We see that 
(anti)hyperons are indeed more abundant than non-strange (anti)baryons. 
Taking into account the disintegration of strange baryons, we are finding a 
much greater number of observed protons dp/dy\°^ ~ 65 ± 5 in the central 
rapidity region. It is important when quoting results from table [| to recall 
that: 

1) we have chosen arbitrarily the overall normalization in table [|, only 
particle ratios were computed, and 

2) the rapidity baryon density relation to rapidity proton density is a 
consequence of the assumed value of X q , which we chose to get E/B ~ 100 GeV 
per participant. 



6.2 Comparison of theory and experiment 

Our study of RHIC results has just begun and we expect that it remains an 
ongoing project as we enter the period of RHIC data collection and analysis. 
Some of our predictions can already be compared with experiment (STAR at 
RHIC). These results are in agreement with our model of how QGP evolves 
at RHIC. 

In order to compare with experiment, we need to adjust from the as- 
sumed 100+lOOAGeV RHIC energy to the RHIC run energy in 2000 which is 
65+65.4 GeV. We consider the A/A ratio. At 8.6+8.6.4 GeV (corresponding 
to 160.4 GeV the experimental result is A/A ~ 0.12 while at 100+100.4 GeV 
we predict A/A = 0.89 ± 0.02 in table ||. We also made a prediction for 
35+35.4 GeV B, where we found A/A ~ 0.5. The interpolated result for the 
65+65.4 GeV RHIC run is A/A ~ 0.68 ± 0.05, which is nearly exactly the 
October 2000 APS-meeting STAR ratio A/A = 0.7 ± 0.05 ± 0.2 . From our 
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A/ A ratio one obtains 

V A^A' 

and this implies taking note of the heading and caption in table ^| that the 
p/p ratio is about 7% smaller compared to A/A ratio. In fact STAR reported 
p/p = 0.65 ± 0.05, in agreement with our result. 

Both particle ratios can thus be well described within our understanding 
of the hadronization and properties of QGP at RHIC. In their absolute values 
reflect on the expected small baryon content at the initial RHIC energy. Both 
should increase by 25-30% at the top RHIC energy (100+100 AGeV). 

The most interesting result seen in table ^, which confirmation we antic- 
ipate and await with anxiety is the hyperon-dominance of the baryon yields 
at RHIC, a fact which as we believe does not depend on a detailed model 
hypothesis. Another interesting property of the hadronizing hot RHIC mat- 
ter, as seen in table [|, is that strangeness yield per participant is expected 
to be 13-23 times greater than seen at present in Pb-Pb interactions at SPS 
energies, where we have 0.75 strange quark pairs per baryon. 

7 Summary and conclusions 

The enhanced strangeness yield observed at 158-200^4 GeV reactions corre- 
sponds according to our study to 0(0.75) ss-pairs of quarks per participant 
baryon, see section ^. Considering the properties of QGP-liquid explored in 
section |^ and the resulting initial gluon temperatures we have shown that 
this exceptionally high yield is achievable in a short time that the collision 
is known to last by in- plasma gluon-fusion reactions, G + G — > s + s, as was 
proposed many years ago. □ 

Beyond strangeness enhancement, multi strange baryons and antibaryons 
have been recognized as being even a more specific source of information 
about the deconfined state of the matter. The analysis of experimental re- 
sults from 158>lGeV Pb-Pb interactions carried out in sections [| (chemical 
abundance analysis) and in section ^ (spectral m± shape) confirms this. We 
find that strange baryons and antibaryons are dominantly produced in sudden 
hadronization of the quark-gluon fireball with chemical and thermal decou- 
pling occurring at the same time. 

From the appearance of first results we recognized that sudden hadroniza- 
tion is required. EJ However, initially the mechanisms leading to direct 
hadronization (hadron production from QGP fireball without equal rescat- 
tering), have not been theoretically understood. Here, we have presented 



tab20002: submitted to World Scientific on February 8, 2008 31 



in section ^ a complete reaction picture consistent with QGP equations of 
state. Il3 The sudden breakup (direct hadronization) into final state particles 
occurs as the QGP fireball super-cools, and in this state encounters a mechan- 
ical instability, which drives the explosive disintegration. Deep supercooling 
requires a first order phase transition, and this in turn implies presence of a 
latent heat and we have estimated its magnitude in section ^. 

There are several experimental facts that lead us to consider sudden QGP 
hadronization. The most recent point has been the results presented in sec- 
tion ^| which show agreement between thermal and chemical freeze-out anal- 
ysis. An important experimental input is the identity of the m± spectra for 
strange baryons and antibaryons. We have been motivated initially to pur- 
sue sudden QGP hadronization by the need to use chemical nonequilibrium 
to describe hadron production. At first we have proposed this for the case 
of strange quarks, and more recently as the experimental data became more 
comprehensive, we also recognized the need to consider light quark chemical 
nonequilibrium — we described details of these developments elsewhere. Ej 

The light quark pair excess is a convenient way to understand hadron 
multiplicity excess, which we have related to an entropy excess, lj But is this 
the only good way to understand the experimental hadron production data? 
It is hard to check all possibilities. However, we explored (in medium) changes 
of hadron masses. We have found that invariably the statistical significance 
of the analysis decreases as we modify individual hadron properties, as long 
as we have been assuming chemical equilibrium of hadron abundances, and 
a consistent change of both mesons and baryons properties. Our extensive 
trial and error searches have convinced us that the only theoretical descrip- 
tion of hadron production data that works (x 2 /dof < 1) requires excess of 
valence light quark pair abundance, irrespective of the detailed strategy of 
data analysis. 

So how it can be that there are many people who believe that one 
can describe hadron production data assuming chemical equilibrium of light 
quarks. In fact, if we take hadronization to somehow occur in chemi- 
cal equilibrium, or/and ignore the possibility of plasma formation, we find 
% 2 /dof > 2.5, dof = 10 in our analysis^ which agrees when carried out with 
this assumption with results of others. EJ Even though this \ 2 IS without con- 
vincing physical significance, what has caused that this result has not been 
forgotten is that the chemical freeze-out temperature for equilibrium condition 
comes out to be near to the expected equilibrium phase transition condition, 
T ~ 170 MeV, see section |[ Thus the assumption of equilibrium seems to 
be consistent with the properties of QGP phase. Furthermore, the study of 
soft pions suggests (unlike it is with hard pions or other hadrons) that the 
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thermal freeze-out occurs at about 120 MeV. Naturally, soft pions can be pro- 
duced by more complex collective processes not considered in the conventional 
thermal analysis and we thus did not introduce these particles into our study 
preferring to depend on the numerous other results in section ||. 

Further study of details of chemical freeze-out, involving consideration of 
strong resonances may lead to a direct experimental differentiation of these 
two alternatives. We are presently considering S*(1385), eil which if freeze- 
out occurs at T = 170 MeV would be the dominant source of A considering 
the decay £*(1385) — > A + ir (88%). But we are already today quite certain 
about the correctness of the sudden hadronization picture since: 

1. the experimental strange and non-strange, hadron abundance results sta- 
tistically strongly favor chemical nonequilibrium in hadronization; 

2. baryons and antibaryons have the same m± spectra; and 

3. our recent finding from the analysis of m± spectra described in section ^| 
that the thermal freeze-out temperature is the same as the chemical non- 
equilibrium freeze-out condition. 

Our enhanced understanding of the fireball supercooling and resulting 
sudden breakup, as well as our study of thermal freeze-out properties has in 
the past year considerably strengthened our case. Not only we can now argue 
that the point of instability agrees with the thermal and chemical freeze-out, 
but also we have considered the freeze-out surface dynamics and have shown 
that the break-up velocity Vf is nearly velocity of light, as would be expected 
in a sudden breakup of a QGP fireball. 

What are the possible causes of the chemical non-equilibrium? Certainly 
even if the breakup of the QGP-fireball is sudden we could have abundance 
equilibrium. In fact, the microscopic processes governing the fireball breakup 
determine how the physical and statistical properties of the fireball change 
at the breakup point. In particular as gluons convert into quark pairs and 
hadrons, the gluon abundance parameter 7 g — > and the quark chemical oc- 
cupancy parameters have to increase significantly, as we determine in data 
analysis. On the other hand the energy E and baryon content b of the fire- 
ball are conserved. Entropy S is conserved when the gluon content of a 
QGP fireball is transformed into quark pairs in the entropy conserving pro- 
cess G + G — ► q + q. Similarly when quarks and antiquarks recombine into 
hadrons, entropy is conserved in the range of parameters of interest here. Thus 
also E/B and S/B is conserved across hadronization condition. The sudden 
hadronization process also maintains the temperature T and baryo-chemical 
potential /it across the phase boundary. 
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We have described strange hadron m^-spectra with astounding preci- 
sion in all available centrality bins. The remarkable similarity of m± spectra 
reported by the WA97 experiment is interpreted by a set of freeze out pa- 
rameters, and we see that production mechanism of A, A, and S, 5 is the 
same. This symmetry, including matter-antimatter production is an impor- 
tant cornerstone of the claim that the strange antibaryon data can only be 
interpreted in terms of direct particle emission from a dcconfmed phase. In 
presence of conventional hadron collision based physics, the production mech- 
anism of antibaryons is quite different from that of baryons and a similarity 
of the m± spectra is not expected. Moreover, even if QGP is formed, but a 
equal phase of confined particles is present, the annihilation of antibaryons 
in the baryon rich medium created at CERN-SPS energy would deplete more 
strongly antibaryon yields, in particular so at small momenta, with the mnre 
abundant baryons remaining less influenced. This effect is not observed. Ej 

Similarity of m_i_-spectra does not at all imply, and one should not expect 
a similarity of particle rapidity spectra. As hyperon are formed at the fireball 
breakup, any remaining longitudinal flow present among fireball constituents 
will be imposed on the product particle, thus A-spectra containing potentially 
two original valence quarks are stretched in y, which A-y-spectra are not, as 
they are made from newly formed particles. All told, one would expect that 
anti-hyperons can appear with a thermal rapidity distribution, but hyperons 
will not. But both have the same thermal-explosive collective flow controlled 
shape of mj_-spectra. 

The only systematic disagreement we see is that the ft + results have 
a noticeable systematic low p±_ enhancement anomaly visible in two lowest 
jnj_ bins in all centrality bins. This result shows that it is not a different 
temperature of freeze-out of fl + fl that leads to more enhanced yield, but a 
soft momentum secondary source which contributes almost equal number of 
soft f2 + f2 compared to the yield expected in view of the behavior of other 
strange hadrons. 

We have presented here a short overview of our recent work developing 
further strangeness as a signature of the formation of quark-gluon plasma, and 
as a diagnostic tool allowing the study of its properties. The understanding 
of the properties of the QGP formed in CERN-SPS experiment has greatly 
improved in past 18 months. 
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